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Abstract -Motivated by the dynamics of cultural change and diversity, we generalize the three- 
species constrained voter model on a complete graph introduced in [J. Phys. A 37, 8479 (2004)]. 
In this opinion dynamics model, a population of size A'' is composed of "leftists" and "rightists" 
that interact with "centrists": a leftist and centrist can both become leftists with rate (1 + q)/2 
or centrists with rate (1 — g)/2 (and similarly for rightists and centrists), where q denotes the 
bias towards extremism (g > 0) or centrism (g < 0). This system admits three absorbing fixed 
points and a "polarization" line along which a frozen mixture of leftists and rightists coexist. In 
the realm of Fokker-Planck equation, and using a mapping onto a population genetics model, we 
compute the fixation probability of ending in every absorbing state and the mean times for these 
events. We therefore show, especially in the limit of weak bias and large population size when 
\q\ ~ A'^^^ and N ^ 1, how fiuctuations alter the mean field predictions: polarization is likely 
when g > 0, but there is always a finite probability to reach a consensus; the opposite happens 
when g < 0. Our findings are corroborated by stochastic simulations. 
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Introduction. — Understanding how diversity is 



maintained and how traits copied by imitation spread are 
T— I central issues in genetics, ecology, and in behavioral sci- 
ence [IHZ]- In this context, there has recently been an 
ijpsurge of interest in statistical physics models predicting 
biological and cultural change, see e.g. [71IH], with relevant 
^""^ phenomena described by closely related models [3] ■ 

^ • One of the basic issues in opinion dynamics is to un- 
derstand the conditions under which consensus or diver- 



sity is reached from an initial population of individuals 



X 

(agents) with different opinions. The voter model [TU] is 
" " arguably the simplest and most popular opinion dynamics 
model [7]. In the voter model and in its variants [8], the 
evolution is implemented by allowing each agent, viewed 
as a "spin" [TU] , to adopt a new state in response to opin- 
ions in a local neighborhood. While the classic 2-state 
voter model unavoidably evolves towards consensus, it has 
recently been proposed that the competing features of con- 
sensus and incompatibility could be two realistic ingredi- 
ents to help explain cultural diversity as an alternative to 
consensus [S1|B|. The basic idea is that agents with sufS- 
cicntly disparate opinions do not interact breaking up into 
distinct cultural states and in this case no consensus can 



be reached ("incompatibility"), while individuals sharing 
close opinions may evolve towards a global consensus. In- 
fluential examples of models characterized by consensus 
and incompatibility are the Axelrod [5] and the bounded 
compromise models [6] that describe the formation and 
evolution of cultural domains. Recently, a (symmetric) 
discrete three-state version of the bounded compromise 
model was studied in low dimensions (where it exhibits 
slow non-universal kinetics) and solved analytically in 
its zero-dimensional formulation |12j . In such a three-state 
opinion formation model, there are two species, A and 
B, respectively called "leftists" and "rightists", that do 
not interact among them (incompatibility) |ll[|12j . How- 
ever, A and B individuals interact with the third species, 
C ("centrists"), and thus indirectly compete to impose a 
consensus. Due to the A and B incompatibility, the final 
state can be either consensus or polarization with a frozen 
composition of leftists and rightists. 

It is natural to generalize the three-state constrained 
voter model of Refs. [HJdl] by assuming that the inter- 
action between "extremists" {A and B individuals) and 
centrists is characterized by a bias q: extremists are more 
persuasive when 9 > 0, while centrists prevail when g < 0. 
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Here, our goal is to study how fluctuations alter the pre- 
dictions of the mean field rate equations concerning the 
system's fate in the presence of a small bias q. We thus 
determine the "fixation probability" jT3] of each absorb- 
ing state of the system (comprising N individuals on a 
complete graph) and compute the average times for these 
events to occur (mean fixation times). In fact, while polar- 
ization is generally the most probable outcome when g > 
and ccntrism is likely to prevail when q < 0, we carefully 
analyze the effects of the bias on the system's fate, espe- 
cially when the population size iV ^ 1 is large and the 
bias \q\ <^ 1 is weak but Nq = 0(1), and determine the 
finite probability that the final state is a consensus when 
q > and the frozen stationary state when q < 0. 

The 3-state constrained voter model. We con- 
sider a population of N individuals on a complete graph, 
Na arc of species A, Nb of type B and Nc of species C, 
with N = Na + Nb + Nc ■ The idealized complete graph 
is used because of its simplicity and analytical tractabil- 
ity. In the language of the voter model, the species A and 
B represent "radical opinions" (e.g. leftists and rightists) 
while the type C stands for an "intermediate state" (e.g. 
centrists). Therefore, A and B individuals interact with 
C but do not interact among them. Hence, species A and 
B both strive to spread at the expense of C. This system 
therefore evolves through the interactions that species C 
has with types A and B. The latter follow the general 
prescriptions of the voter model and proceed by imitation: 
one individual is picked randomly and adopts the opinion 
of one of its random neighbor provided that at least one 
of the individual is of species C. In this way, the system's 
dynamics at each time increment can be schematically de- 
scribed by the following reactions: 
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where — 1 < g < 1. The special case q = was thoroughly 
studied in Ref. [12] and will therefore not been discussed 
here. The parameter q measures the bias towards polar- 
ization (when g > 0) or centrism (when q < 0). In the 
former situation extremists {A's and i?'s) are more per- 
suasive than centrists (C"s), while in the latter centrism is 
the dominating and more persuasive opinion. As a conse- 
quence, we anticipate that polarization is the most prob- 
able final state when g > 0, while centrism consensus is 
expected to be the most likely stationary state when q < 0. 

At mean field level, assuming a population of infinite 
size {N — ^ oo) and the absence of any random fluctuations, 
the system dynamics is described by the following rate 
equations (REs) for the densities a = Na/N and b = 



Nb/N of species A and B, respectively: 



dt 



a = qa{l — a — h), —b = qb{l — a — b), 



(2) 



where we have used the fact that a{t) + b{t) + c{t) = 1. 
The REs ^ admit 3 absorbing fixed points, (a, 6, c) = 
(a(oo),5(oo),c(oo)) e {A ^ (1,0,0), 6 (0, 1,0),C = 
(0,0, 1)} and a line of fixed points given by AB = (a, 1 — 
a, 0), with < a < 1. In fact, the system ^ can be solved 
exactly, yielding [H] 



a{t) 



=9* 



1- {x + y){l-ei*) 



b{t) = - 



ye 



qt 



[x + y){l - e'i^)'' 



where x, y respectively are the initial densities of species 
A and B, i.e. x — a(0) and y = 6(0). These results reveal 
that the line of fixed points is the stable solution when 
g > 0, i.e. in this case (a, b, c) = (xj [x + y), yjix y), 0), 
with a -I- 6 = 1. On the other hand, when g < and the 
bias favors the species C, one finds a = 6 = and c = 1. 
Furthermore, we notice that the REs ^ predict that the 
ratio of the densities is conserved, i.e. a{t)/b{t) = x/y. 

When the population size N is finite, demographic fluc- 
tuations can drastically alter the mean field predictions. 
In such a setting, the dynamics is no longer determin- 
istic and the system's fate depends non-trivially on the 
bias strength and on the population initial composition, 
parametrized by the initial densities {x,y). Within a 
stochastic formulation of the model, the moves ([Ij define a 
birth-death process [14], where, according to the num- 
ber of individuals of each species increases or decreases by 
one unit in each time step. A quantity that is central 
for our discussion is P'^'^{x,y), the polarization fixation 
probability along the absorbing line AB starting from an 
initial population composition (x,y,l — x — y). Hence, 
pAB gjygg tj^g probability to find the system locked into 
a polarized state where "extremists" {A's and B's) coexist 
without interacting. This probability obeys the following 
backward master equation [14j : 



(r+ + T- + T+ + T-)P-^«(x, y) - T~P 



T+P-^'^ix + <5, y) + r-p-^«(x, y-S)+ T+P-^«(.T, y 



where, in the limit iV ^ 1, the transition rates are 
r± = (l±g)x(l-x-2;)/2andT± = {l±q)y{l-x-y)/2, 
with S = N~^. This two-dimensional equation has to be 
supplemented by the boundary conditions: P'^^{x,0) = 
P^B(0,y) = and P^^(x, 1 - x) = 1. By Taylor- 
expanding ^ to second-order in S, one finds 



s[xda: + ydy] + 



xdl 



jAB 



{x,y)^o, (4) 



where we have introduced the parameter s = Nq. When 
\s\ ^ 1, the diffusion term on the left-hand-side of (H]) is 
negligible in front of the deterministic drift term, while the 



p-2 



opposite occurs when \s\ <^ 1. This implies that the in- 
teresting situation arises when s is of order one, i.e. when 
N ^ 1 and \q\ ~ N^^; otherwise one would essentially 
recover the mean field predictions of ([2|) when \s\ ^ 1, or 
the results of Ref. [TJ when |s| <C 1. It is also useful to 
notice that the associated backward Fokkcr-Planck (FP) 
differential operator £bFP reads [H] : 



C 



bFP 



(1 



y) 



2N 



2s{xd^ + ydy) + xdl + ydl]{5) 



A one- dimensional version of £bFP often appears in pop- 
ulation genetics [1] and evolutionary game theory [5| . 

Fixation Probabilities. — We have seen that the 
mean field treatment ([2]) predicts that the line AB and the 
fixed point C are the system's attractor when q > and 
q < 0, respectively. The quantity P-^'^ , that obeys 
and the probability density F^^{x,y) that the system's 
final state has coordinate (a, 1 — a) along the absorb- 
ing line, are central to study the system's fate in the 
presence of fiuctuations. These quantities are related by 
P^'^{x,y) ~ daF^'^{x,y). Hence, F^'^{x,y) obeys 
the same backward FP equation ^ as p-^^ but with 
the boundary conditions P^^(0,?/) = F^^{x,0) = and 
F^'^{x, 1-x) = 5{a-x). As in Ref. HI], Eq. g]) turns out 
to be separable and exactly solvable. To obtain its solution, 
it is useful to introduce the polar coordinates (p, 9) such 
that y/x = pcos9 and ^/y ~ psm9. In these coordinates, 
with < p < 1 and < 6* < 7r/4, Eq. (g]) becomes: 



[{4sp-p-')d, + dl] P-^''ip,9) 
+ p-^{{ta.n9-cot9}de + d^)P-^'^{p,9) 



(6) 



0. 



with boundary conditions P^'^ [p = Q,9) = and 
P-^>^{p = 1,6*) ^ 1. The probability density F^^ also 
obeys the equation ([5]) but with the boundary conditions 
F^^{P = 0,6*) = and F^'^{p ^ 1,9) = 5{a - cos^9). 
Following Ref. [12], we write P^'^ and F^'^ in the form 
X^n CnRn{p)un{9) sin (20), whcrc the Rn{p) and Un{9) are 
the eigenvectors (with eigenvalues A„ ) of the following ra- 
dial and angular Sturm-Liouville problems, respectively: 
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(sin^ 



[Asp^ - 1] - X„Rn = 



(1 + Xn)Un 



(7) 







The solution to the radial equation ([7]) with 
boundary condition i?„(0) = reads Rn{p) c(. 



sp p 2ri i/„^^y2(sp^), where /„ is the modified Bessel 



function of first kind and order n [TS]. The equation 
([7]) for the angular function u{9) coincides with a sta- 
tionary Schrodinger equation in a Poschl- Teller potential 
hole whose solution is u„(0) cx ^sin (26') J'^^;^ (cos (26')), 
where P^ denotes an associate Legendre polynomial of 
first order [T^[Tn]. The eigenvalues are found to be 
A„ = 4(n -I- l)(n -|- 2) and the coefficients c„ are deter- 
mined using the orthogonality of the P^'s together with 
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Fig. 1: (Color online). Fixation probabilities for s > as 
functions of a;: (o), (A), (□); P-^^ (o) is compared 
with © (solid curve). Parameters are TV = 200, s = 4 (i.e. 
q = 0.02). (a) Initial density of yl is a; = y and P^ ^ P'^ . 
Inset: Comparison of P^ with its analytical approximation, 
see text, (b) Initial density of y4 is a: = 2y. Inset: final densities 
of species A {x) and B (-(-) as functions of a; = 2y. Numerical 
results have been averaged over 2 x 10^ samples. 



the boundary conditions for F^'^ and P^^ 
leads to 



This 



In+l/2{s{x + y)) 



j£y_s{i-x-v) 
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2n+ 1 



-1/2 
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^ \n{n+l) 
x-y\ P,K2a-I) 



(8) 



jx + yj ^a{l - a) 

Since P"^^ ix,y) = da (x, y) , using the properties of 
the associate Legendre polynomials [TS], one also obtains 



P-^S(x,y) = 2 



Xy_ s(l-x~y) 



y 
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n odd 



2n + 1 
nin + 1) 



In+i/2{s[x + y)) y pi ( x~y 



1+1/2 



x + y 



(9) 



where the subscript "n odd" means that the sum runs over 
odd integers n. Since p-^^ is the polarization fixation 
probability, the probability that the system's final state is 
consensus (with either A, B, or C) is 1 - P^'^ . 

The quantities P-^^ and the fixation probabilities of the 
absorbing states A, B and C, respectively denoted P^, P^ 
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and P'^, are related by P^{x,y) + P'^ {x,y) + P^'^ {x,y) = 
a+b = 1—P'^{x, y). In fact, as C's interact identically with 
A's and S's, a mapping onto a population genetics model 
where extremists are regarded as mutants of a single class 
yields P'^ = - — jj^"!^/ — (see below). Furthermore, the 
fixation probabilities and the species density are related by 
by a = P^{x,y) + da' a' Fa- {x , y) and b = P'^{x,y) + 
P^^{x,y)- J^da'a'Fa-{x,y), with P'^{x,y) ^ P-^iy,x). 
When \s\ < 1, /„+i/2(s(a; + 2/))//„+i/2(s) « + 



1.0' 



y) 



i+l/2 



[T5] and, 



with ([5]) and (O, one obtains 

Loi^.y) - Pm^) ^ " '^"^ + y - 1} + From 
this expression and the compact result P;^q{x,x) = 1 — 
-^^T obtained in [ig, we infer P^^{x,x) « 6(1 + s)x'^ 
when X = y 1. Hence, at low initial density x = y and 
for weak bias ( | s | <C 1 ) , P'^^ is a quadratic polynomial in 
X with an amplitude proportional to s. 

Fixation probabilities when s > 0. When s > 0, we 
find that the fixation probability P"^^(a;, y) displays a sig- 
moid shape, that steepens when s increases, interpolating 
monotonically between (when x -\- y ^ 1) and 1 (for 
x + y — >■ 1), see Fig.[TJ The agreement between the analyt- 
ical prediction © and the results of stochastic simulations 
(obtained using the Gillespie algorithm [T^) is excellent, 
as illustrated in Fig. [1] This figure also shows that P'' 
steeply decays to zero when x + y increases, as expected 
from its analytical expression. Furthermore, P^ and P^ 
display positive skewness and maxima around the initial 
densities (a;* , y* ) such that P-^^ (x* , ) = P'^ (.t* , ) 
(x* + y» w 0.157 in Fig. [Ha)). Moreover, as shown in 
the inset of Fig. [Ijb), the stationary densities steadily 
approach their mean field values, i.e. a ^ xjix: + y) and 
b — > x/{x + y), when a; + y > 0.30. In such a regime, 
one thus has a w ^ = P^{x,y) - da' a' Fa' {x , y) 
and, with ([5]), this yields the following approximation 
for the fixation probability of A: P-^{x,y) « — 

V^e {n{n+l) \ /„+i/2(s) / ) ^^'^ 

where A = {x — y)/{x + y). The inset of Fig. [Tfa), 
shows that the agreement between this approxima- 
tion and the results of numerical simulations increases 
with X = y > 0.20. A similar approximation gives 
pBr:^l_ P^^{x, y) - P^{x, y). 

Fixation probabilities when s < 0. When s < 0, the 
bias is towards the absorbing state C = (0, 0, 1) and the 
fixation probabilities P^ and P^ are vanishingly small. 
For small initial densities of A and B (i.e. x-l-y <C 1), the 
probability of reaching the absorbing line AB is also very 
small (< N~^). However, this fixation probability P^'^ 
grows monotonically when x + y increases and, according 
to ([9|) and as shown in Fig. [2ja,b), approaches the value 
one when x-\-y — >■ 1. Figure[5]shows that the system's most 
likely fate is either to end up in the state C or on the line 
AB. To determine a concise and accurate approximation 
of p-^'^ , one can use a mapping onto a population genetics 
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Fig. 2: (Color online). Same as in Fig. [1] with s = —4 in 
semi-log scale, and comparison with ((9| (solid curve) and (|10|l 
(dashed curve), see text. Tlie initial density of ^ is a:: = y in 
(a) and a; = 2y in (b). Inset of (a): P'' and P'^'^ in linear 
scale. Insets of (b): stationary density of ^ (x) as function of 
X — 2y compared with the analytical expression (solid curves) 
given in the text, in linear (left) and semi-log (right) scales. 



model where C is considered as a "wild-type" allele and 
both A and B species are regarded as forming a single- 
class deleterious allele [T] . In this situation, the probability 
P-^'^{x, y) = P-^^(a;-|-y) corresponds to the fixation of the 
(single-class) "mutants" and depends only on a; -I- y and 
s. The angular dependence therefore drops out from ^ 
and, with the (new) boundary conditions P-^^(a; + y = 
0) = and P^^{x -t- y = 1) = 1, one obtains the following 
approximation for P^'^ ; 



P-^^{x + y) 



(10) 



e2k| _ 1 ■ 

As shown in Fig. [21 this result is found to be an excellent 
approximation of P^^ when a:+y > 0.30. In such a regime 
P-4 ss pB « and P^^{x + y) I - P'^ {x + y) (see inset 
of Fig. [H^a)). Using (|TU|) . one can obtain the stationary 

densities b ~ a(y/a;) and a ~ (^) (^^^^m^rr^)- The 
excellent agreement between this expression and numerical 
simulations is demonstrated in the insets of Fig. [5] (b) . 

Mean Fixation times. — The mean times necessary 
to reach the absorbing states, also called the mean fixation 
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times (MFTs), are quantities of great interest in evolution- 
ary dynamics [I1[1I11[71|S] . Here, the unconditional MFT 
to reach any of the system's absorbing states is a function 
r(x, y) of the initial density (x, y) obeying the backward 
FP equation C}:iYY'{x,y)T{x,y) — —1 [see ^] with bound- 
ary conditions t(1, 0) = r(0, 1) = r(0, 0) = r(a, 1 — a) = 
[Til 14], Such an equation can be solved by using an ex- 
act mapping onto a suitable population genetics model. 
In fact, since we are interested in the unconditional MFT, 
and species C docs not make any distinction between A 
and B individuals, as seen above, the latter can be con- 
sidered as belonging to a single "mutant type" class that 
interacts with the "wild type" C [T]. In this formulation, 
the equation for r becomes one-dimensional with two ab- 
sorbing states: either only wild-type individuals (consen- 
sus with C) or only mutants (either A ov B consensus, 
or polarization in a frozen mixture of A's and i?'s). The 
unconditional MFT thus depends only on the initial den- 
sity a; -I- y of extremists (A's and S's): r(x, y) = t{x + y). 
Hence, with the variable w = x + one obtains: 



N 



2s 



dr^w) d^T{w) 
dw dw'^ 



-1, 



(11) 



with t(0) = r(l) = 0. Eq. pT|) frequently appears in pop- 
ulation genetics, where it describes the MFT in a diallelic 
haploid population in the presence of selection of (weak) 
intensity \q\ = \s\/N ^ 1 [1] (see also [4]). An interesting 
property of is its invariance under the transformation 
(w, s) — > (1 — w, — s). This implies that t{w) for s > co- 
incides with r(l — w) for s < 0, see Fig. [31 Equation ([TT|) 
can be solved by standard means (see e.g. [11[I1] and ref- 
erences therein) yielding a cumbersome expression. The 
latter takes a more compact form for the special value 
w = 1/2, when it reads: 



N 



s(l + e*) 
N 



e^''Ei(-2s) + (e'* - l)ln(2|s|; 
[e^(e^-l)7E + Ei(.s)-Ei(2s)] 



2se'*(l + e'') 

+ 2g(l + e-) [2Ei(^)-2e-(e^ + l)Ei(-^)]. (12) 

where Ei(a::) = jHoo T denotes the exponential integral 
and 7e = 0.5772... is Euler-Mascheroni's constant [T5] . 
From the general expression of r(a; -|- y), one finds that 
the unconditional MFT scales linearly with the system 
size N, i.e. t{x + y) ~ Nfr{s,x + y), where fr is a 
scaling function. Such a scaling relationship is fully con- 
firmed by the numerical results reported in Fig. [2Ka,b). 
We can infer from the latter that fr has an inverted u- 
shape dependence on w ^ x + y and skewness towards 
small values of w for s > and large values of w for 
s < 0, see Fig. [2Ka,b). This is a consequence of the 
{w, s) — ?> (1 — w, —s) symmetry obeyed by t{x + y). The 
conditional mean fixation times, r"^, to reach the absorb- 
ing states S G {A,B,C,AB) can be obtained by solving 
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Fig. 3: (Color online). The normalized unconditional mean 
fixation time (MFT) t/N (x), compared with the solution of 
Eq. (solid curve) , and the (normalized) conditional MFTs 
T-^/N (o), r^/iV (•), T^/N (□), T^^/N (o), see text. The 
parameters are A'^ = 200, x — y (initial A density), s = 4 in (a) 
and s = —4 in (b). The maxima of the unconditional MFTs 
are r(0.22) ^ OMN in (a) and r(0.78) « 0.94iV in (b). The 
numerical results have been averaged over 2 x 10^ samples. 



Ci,Fp{x,y)[P^{x,y)T^{x,y)] = ^P^ {x,y), where P^ix,y) 
is the fixation probability of the state S, with the appro- 
priate boundary conditions [I1[I1]. Except for r'', these 
equations cannot be mapped onto one-dimensional prob- 
lems and are difficult to solve. However, as p-^'^ — > 1 when 
w —> 1 and s > 0, we infer that in this regime the equa- 
tion for T-^'^ coincides with ([TT|) and therefore t-^'^ ~ t, 
as confirmed by Fig. [3l^a). Similarly, since P'' -> 1 when 
i(j <C 1 and s < 0, coincides with r in such a regime, see 
Fig. [3^b). The numerical results reported in Fig. [31Ia,b) 
show that all r'^'s scale linearly in N (and monotonically 
with a; + y), in a manner that appears to be independent 
of the sign of s. In fact, while the results for and in 
Fig-Ob), where s < 0, and those of Fig.[3la), where s > 0, 
are comparable, the former are more noisy than the lat- 
ter. Furthermore, as also found for the special case q = 
[T^ . the extremists MFTs, i.e. and t^, are always the 
longest mean fixation times. 
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Conclusion. — Wc have generalized a basic three- 
state opinion dynamics model introduced in Refs. pTlfT^ 
and considered a finite population of N individuals con- 
sisting of leftists, rightists and centrists interacting on a 
complete graph. Motivated by recent studies concerning 
the formation of cultural diversity, the system's dynamics 
is characterized by two competing features: (i) "extrem- 
ists" (leftists ^'s or rightists B's) interact with centrists 
(C's) and at each elemental step an extremist can become 
a centrist with rate (1 — '7)/2, while a centrist can become 
either leftist or rightist with rate (1 -I- q)/2; (ii) extrem- 
ists do not interact. The former feature drives the system 
towards consensus while the latter accounts for incompat- 
ibility and leads to a frozen steady state of A's and B's 
(polarization). The parameter q denotes a bias favoring 
polarization when g > and centrism when q < 0. 

This three-species voter model is characterized by three 
absorbing fixed points (one for each species) and an ab- 
sorbing line where a frozen mixture of extremists coexist 
(polarization). Here, we have studied the model's final 
state properties and showed that fluctuations drastically 
alter the mean field predictions in the presence of a small 
bias. In fact, while polarization is generally the most prob- 
able outcome when q > there still is a finite probability 
to attain a consensus, and the opposite situation arises 
when q < 0. Our results are particularly relevant in the 
limit of large population size N ^ 1 and vanishing bias 
l^l ^ 1 with s = Nq = C(l) finite, where the fluctuations 
and the deterministic drift are of the same intensity. Such 
a situation corresponds to the "weak selection limit" fre- 
quently considered in life and behavioral sciences [Tl|3l|4]. 

The polarization fixation probability, P-^'^ , has been 
obtained analytically by solving a (separable) backward 
Fokker-Planck equation and our results have been checked 
against stochastic simulations. Via a mapping onto a pop- 
ulation genetics model, we have shown that when s < 
P"^^ approximately grows exponentially with the initial 
density x + y of extremists, whereas the centrist fixation 
probability P'' decays (approximately) exponentially with 
X + y when s > 0. The mean fixation times (MFTs) have 
been computed and found to scale linearly with N (see 
Fig. [3]). Furthermore, the unconditional MFT r has been 
shown to satisfy t{x + y,s) = t(1 — x — y, —s). 

Our results show that mean field rate equations cannot 
describe the final state of a simple opinion-dynamics model 
in the presence of a small bias. This further illustrates the 
pertinence of statistical physics methods to describe the 
evolutionary dynamics of models of cultural diversity. 
Acknowledgments: Sid Redncr is gratefully acknowl- 
edged for a useful discussion. 
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